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ABSTRACT 

This work presents blind joint interference suppression and power 
allocation algorithms for DS-CDMA networks with multiple relays 
and decode and forward protocols. A scheme for joint allocation 
of power levels across the relays subject to group-based power con- 
straints and the design of linear receivers for interference suppres- 
sion is proposed. A code-constrained constant modulus (CCM) de- 
sign for the receive filters and the power allocation vectors is devised 
along with a blind channel estimator. In order to solve the proposed 
optimization efficiently, an alternating optimization strategy is pre- 
sented with recursive least squares (RLS)-type algorithms for esti- 
mating the parameters of the receiver, the power allocation and the 
channels. Simulations show that the proposed algorithms obtain sig- 
nificant gains in capacity and performance over existing schemes. 

Index Terms — DS-CDMA, cooperative systems, optimization 
methods, blind algorithms, resource allocation. 

1. INTRODUCTION 

Multi-antenna wireless communication systems can exploit the spa- 
tial diversity in wireless channels, mitigating the effects of fading 
and enhancing their performance. Due to size and cost of mobile 
terminals, it is usually impractical to equip them with multiple anten- 
nas. However, spatial diversity gains can be obtained when single- 
antenna terminals establish a distributed antenna array via coopera- 
tion (i|-|2|. In a cooperative system, terminals or users relay signals 
to each other in order to propagate redundant copies of the same sig- 
nals to the destination user or terminal. To this end, the designer 
must use a cooperation protocol such as amplify-and-forward (AF) 
t2l and decode-and-forward (DF) l2ll3l. 

The use of cooperative diversity and multiple hops is key for 
networks that need to increase the link reliability and extend their 
coverage |2|. Prior work on cooperative multiuser DS-CDMA net- 
works has focused on the assessment of the impact of multiple access 
interference (MAI) and intersymbol interference (ISI), the problem 
of partner selection 13] |6], the bit error rate (BER) and outage per- 
formance analyses f7), resource allocation (4, 5 1 and training-based 
joint power allocation and interference mitigation strategies 1 1 1 , 12) . 
However, these strategies require a significant amount of training 
data and signalling, decreasing substantially the spectral efficiency 
of cooperative networks. This problem is central to ad-hoc and sen- 
sor networks 1 13 1 that employ spread spectrum systems and multiple 
hops. This calls for methods to decrease the amount of signalling and 
training in cooperative wireless networks. 

In this work, blind joint interference suppression and power allo- 
cation algorithms for DS-CDMA networks with multiple relays and 
DF protocols are proposed. A blind scheme that jointly considers the 
power allocation across the relays subject to group-based power con- 
straints and the design of linear receivers for interference suppres- 
sion is proposed. The idea of a group-based power allocation con- 
straint is shown to yield close to optimal performance, while keeping 
the signalling and complexity requirements low. A code-constrained 
constant modulus (CCM) design 1 16|-|20| for the receive filters and 
the power allocation vectors is developed along with a blind channel 
estimator. The CCM design is adopted as it achieves a performance 
close to training-based algorithms. In order to solve the proposed 



optimization problem efficiently, an alternating optimization strat- 
egy is presented with recursive least squares (RLS)-type algorithms 
for estimating the parameters of the receiver, the power allocation 
and the channels. 

The paper is organized as follows. Section 2 describes a co- 
operative DS-CDMA system model with multiple relays. Section 3 
formulates the problem, the CCM design of the receive filters and the 
power allocation vectors subject to a group-based power allocation 
constraint, and a blind channel estimator. Section 4 presents the al- 
ternating optimization strategy along with RLS-type algorithms for 
estimating the parameters of the receiver, the power allocation and 
the channels. Section 5 presents and discusses the simulation results 
and Section 6 draws the conclusions of this work. 

2. COOPERATIVE DS-CDMA NETWORK MODEL 

Consider a synchronous DS-CDMA network with multipath chan- 
nels, QPSK modulation, K users, A'^ chips per symbol and L as the 
maximum number of propagation paths for each link. The network 
is equipped with a DF protocol that allows communication in mul- 
tiple hops using rir fixed relays in a repetitive fashion. We assume 
that the source node or terminal transmits data organized in packets 
with P symbols, the system can coordinate cooperative transmis- 
sions, and the linear receivers at the relay and destination terminals 
are synchronized with their desired signals. The received signals are 
filtered by a matched filter, sampled at chip rate and organized into 
AI X 1 vectors rsd[m.i], Vsn [mi] and rr-dlrrij], which describe the 
signal received from the source to the destination, the source to the 
relays, and the relays to the destination, respectively, 

K 

rsd[mi] = y^ a'^d[mi]Ckhsd,k [mi]bk [mi] 



+ ilsd[mi]+nsd[mi], 

K 

rsTj[mi] = ''^al^^[mi]CkhsT^.k[mi\hk[mi\ 
ft=i 
+ n^rj"?.i]+nsrj[mi], 

K 

•rrjd\mj] = y^al.d\mj]Ckhr-d,k\m.j]hk\mj\ 
ft=i 
+ 'nTjd[mj] +nr^d\mj], 



(1) 



whereM = N + L-l,mj = {j -1)P+1, . . . ,jP,i = 1,...,P, 
j = 1 . . . , rip, P is the number of symbols in the packet, rip — 
n,. + 1 is the number of transmission phases or hops, rir is the 
number of relays, rrij is the index of original and relayed signals, 
nsd[m.i], risr -[m.!] and rir d[mj] are zero mean complex Gaussian 
vectors with variance a^ generated at the receivers of the destina- 
tion and the relays from different links, and the vectors ri^j^lmi], 
rj^^ _ [nil] and 'n^.^[mj] represent the intersymbol interference (ISI). 

The quantities bk[m.i] and bk[mj] represent the original and recon- 
structed symbols by the DF protocol at the relays, respectively. The 
amplitudes of the source to destination, source to relay and relay to 



destination links for user k are denoted by a^^[mi], a'^r- [mi] and 
aj:.jj[mj], respectively. The M x L matrix Ck contains versions of 
the signature sequences of each user shifted down by one position at 
each column as described by 



Cfc(l) 

Cfc(iV) 








(2) 



where c^ = [cft(l), Cfc(2), ..., Cfc(7V)] stands for the signa- 
ture sequence of user k, the L x 1 channel vectors from source to 
destination, source to relay, and relay to destination are Aiad,fc['Tii], 
hsr.-,k[m.i], hrd.klmj], respectively. By collecting the data vec- 
tors in (TJ (including the links from relays to the destination) into a 
{ur + 1)M X 1 received vector at the destination we obtain 



rsd[mi\ 
rrid[m2] 



fr^^dlirin 



''Y^^i a-sd[mi]C khsd,k['mi]bk[mi] 

J2k^l 0-ridlm2]Ckhrid,k[m2]bl.^ [7712] 

-EfeLi a-r„^d['mup]Ckhr„^d,kl'mnp]bl"'' [m„j, 

+ 'n[i] +n[i] 

(3) 



Rewriting the above signals in a compact form yields 

K 

r[i] =J2^i'[i]Ak[i\Ckhk[i]+'n[i]+n[i\ 



K 



(4) 



Y^ Pk[i]Bk[i\ak[i\ + v[i\ + ^tW, 



the links, the {rir + 1) x (n,. + 1) diagonal matrix Ak[i] is given by 
Ak[i] = diag{afe[mi]}, and the (jir + 1)M x (n,. + l)Af diagonal 



power constraints along with a blind channel estimator. To this end, 
the {rir -I- 1)M X 1 received vector in (|4) can be expressed as 

r\i] = Ps\i]Bs\i]as.k[i\ + ^ Pk\i]Bk[i\ak[i\ + v[i] + n[i\, 

k^S 

(5) 
where S = {S\ ,82, ■ ■ ■ , So} denotes the group of G users to con- 
sider in the design. The {rir + 1)M x Gijir + 1) matrix Ps = 
[Psi Ps2 ■ ■ ■ Psq] contains the G effective signatures of the 
group of users. The G(nr + 1) x Gijir + 1) diagonal matrix Bk [i] = 
diag(fes,H yi^[i\...%-''\i] ... 65, W b2^\i]...%l'm con- 
tains the symbols transmitted from the sources to the destination and 
from the relays to the destination of the G users in the group on the 
main diagonal, the G{7ir + 1) x 1 power allocation vector as,k[i] = 



where the (rir + 1)M x (rir + 1)L matrix Ck ~ diag{Cfc . . . Ck} 
contains copies of Ck shifted down by AI positions for each 
group of L columns and zeros elsewhere. The (n^ + 1)L x 1 
vector hk[i] contains the channel gains of the links between the 

source, the relays and the destination, and Pf.[i] = Ckhk[i] is 
the effective signature for user k. The (n,- + 1) x {rir + 1) di- 
agonal matrix Bk[i] = diag(6fe[mi] bl^''[77i2] . . .b'l"'' [rrinj) 
contains the symbols transmitted from the source to the desti- 
nation (bk[i]) and the n,. symbols transmitted from the relays 

to the destination (b^^ [7712] . . . fej.""" [jti„^]) on the main diago- 
nal, and the {rir + 1)M x {rir + l)i\f diagonal matrix Bk[i] = 
diag(bfe[mi] (g) Im bl^''[m2] ® Im ■ ■ ■ bl"'' [m„J (g) Im), where 
(g) denotes the Kronecker product and Im is an identity matrix 
with dimension AI. The {ur + 1) x 1 power allocation vector 

r d [m-np]]'^ has the amplitudes of 



.■Sir 



VidL 



.■Si 



Vfdl 



of 



the amplitudes of the links used by the G users in the group. 



3.1. Blind CCM Receiver Design and Power Allocation Sclieme 
witli Group-Based Constraints 

The linear interference suppression for user k is performed by the 
receive filter Wk[i] = ["'fc.iW' • ■ • > ^"fc,(n,+i)A/M] with {ur + 
1)A4 coefficients on the received data vector r[i] and yields 



ZkW 



Wk[i]r[i] 



(6) 



where Zk [i] is an estimate of the symbols, which are processed by a 

slicei Q{-) that performs detection and obtains bk[i] ~ Q{zk[i]). 

Let us now detail the CCM-based design of the receivers for user 
k represented by Wk [i] and for the computation of the G{nr + 1) x 1 
power allocation vector as.k[i]- This problem can be cast as 



\wT^ , a°J'l] = ai-g min E[{\wk [i]r[i]\' ~ 1)^] 



^k ' "-S.fcJ 



(7) 



subject to as^k[i]0'S,k[i] = Pg and Wk [i]Pk[i] = v, 



where 1/ is a parameter used to enforce convexity |18|. The CCM 
expressions for the receive filter lUk [i] and the power allocation vec- 
tor as,k[i\ can be obtained with the method of Lagrange multipliers 
which transforms (|7) into the Lagrangian function 

Ck = E[(\w'i^[{][Ps[i\Bs[i]as,k[A 



■^Y^Pk[i]Bk\i]ak\i]+r)\i]+n[i])\^ -Vj ] (8) 

k^S 

+ Afc(af_fc[j]a5,fc[i] ~ Pa) + Pk{wk[i]Pk[i] ~ ^), 

where A^ and pk are Lagrange multipliers. An expression for 
0'S,k[i] is obtained by fixing tOfc[i], taking the gradient terms of the 
Lagrangian and equating them to zero which yields 



as,k[i] = (■Rs,fcW + Afc/) ^ds,k[i] 



(9) 



matrix Ak[i] = [a*d[mi] (g) Im flridin^a] (g) Im ■ ■ ■ 
The (ur + 1)M X {rir + 1) matrix Pk has copies of the effective 
signature Pf. [i] shifted down by M positions for each column and 
zeros elsewhere. The {ur + 1)M x 1 vector ri[i] represents the ISI 
terms and the {rir + l)Af x 1 vector n[i] has the noise components. 



where Rs,kli\ = E[\zk[i\fB§[{\P§[i\wk[i\w^[i\Psli\Bs[i]] 
dfwirip] (g Im^ a.G(n,. + l) x G(7ir + 1) correlation matrix and the G(nr + 1) x 1 
vector d5^fc[i] = E[zk[i]B'^[i]Pg[i]wk[i]] is a cross-correlation 
vector. The Lagrange multiplier A^ plays the role of a regularization 
term and has to be determined numerically due to the difficulty of 
evaluating its expression. Now fixing as^k[i], taking the gradient 
terms of the Lagrangian and equating them to zero leads to 



3. PROPOSED BLIND RECEIVER DESIGN, POWER 
ALLOCATION AND CHANNEL ESTIMATION 

In this section, a joint blind receiver design and power allocation 
strategy is proposed using the CCM approach and group-based 



Wk[7] ^ Rk'\i]{dkl^-Pk[^lk'mPk[mk'\i]dkl^-'^), (10) 



where 7fc[i] = p^[i\R^^\i]p^\i] 
by Rk[i\ = E[\zk[i\fr[i\r"[i\] 



, the correlation matrix is given 

and dk\i] = Sfzfcfilrfill is a 



(n,. + 1)M X 1 cross-correlation vector. The quantities Rk[i] and 
dk[i] depend on tlie power allocation vector as,k[i]- The expres- 
sions in ^ and JlOt do not have a closed-form solution as they arise 
from a higher-order optimization. Moreover, the expressions also 
depend on each other and require the estimation of the channel vec- 
tor hk[i]. Thus, it is necessary to iterate ( 110) and (|9) with initial 
values to obtain a solution and to estimate the channel. The network 
has to convey the information from the group of users necessary to 
compute the group-based power allocation including the filter Wk[i]- 
The expressions in ( 110) and (|9) require matrix inversions with cubic 
complexity ( ©(((n^ + l)Mf) and 0{{K{nr + 1))^). 

3.2. Blind Cooperative Cliannel Estimation 

In order to blindly estimate the channel in the cooperative system 
under study, let us consider the covariance matrix R — E[r[i]r^ [i]] 
and the transmitted signal Xk[i] ~ Ak[i]Bk[i]Pk[i]- Let us now 
perform an eigen-decomposition on R 



R 



J2Elxkmx^m + E[ri[^rj"[^\ 



1 2 r 
+ cr I 



(11) 



\^s4>n 



a' I 







where (p^ and (/),j are the signal and noise subspaces, respectively. 
Since 4)^ and <j)^ are orthogonal, we have the condition <p^Xk [i] = 
(f>^Ak[i\Bk[i\Pk[A = (j}"Ak[i\Bk[i\Ckhk[i\ = and hence 



r = /if W Cf Bf [i]Af [i]</.„</.^AfcHBfc[i]Cfc hk[i] 



(12) 



The above relation allows to blindly estimate the channel hk [i] . To 
this end, we need to compute the eigenvector corresponding to the 
smallest eigenvalue of Y^. It turns out that we can use the fact that 
limp_>cx)(Jt/cr^)~'' ~ 4'n'Pn |l^l and, in practice, it suffices to use 
p = 1 or 2. Therefore, to blindly estimate the channel of user k in 
the cooperative system we need to solve the optimization problem 

hk[i] = argmin/if [i]Tfc/ifc[i], subject to ||/ifc[i]|| = 1, (13) 

Alt [i] 

In what follows, computationally efficient algorithms based on an 
alternating optimization strategy will be detailed. 



4. PROPOSED ADAPTIVE ALGORITHMS 

In this section, we develop joint adaptive RLS-type algorithms using 
an alternating optimization strategy for efficiently estimating the pa- 
rameters of the receive filters, the power allocation vectors and the 
channels. Note that the proposed algorithms did not have problems 
with local minima and converge to the desired solutions. 

The first task in the proposed scheme is to build the group of 
G users that will be used for the power allocation and receive fil- 
ter design. A RAKE receiver is employed to obtain z^ ^^H = 



{Ckhk[i\)"r[. 



• Pf. [i]r[i] and the group is formed according to 



compute the G largest \zk [i] 



1,2,..., K. (14) 



The design of the RAKE and the other tasks require channel estima- 
tion. In order to solve ( 113) efficiently, a variant of the power method 
IT6I that uses a simple shift is adopted 



hk[i\ = {I - Tk[i\rk\i])hk[i - I], 



(15) 



where Tk\i] = l/tr[Tk[i]] and hk[i] <- hk\i]/\\hk\i]\\ to nor- 
malize the channel. The quantity Tfe[i] is estimated by Tfe[i] = 



ar[i - 1] + Cf Bfc [{\Ak [i]R ^[i]Ak\i]Bk[i\Ck, where a is a 

forgetting factor that should be close to 1 and R [i] is computed 
with the matrix inversion lemma. The power allocation and receive 
filter design problems outlined in (|7) are solved by replacing the 
expected values in ^ and l llOt with time averages, and RLS-type 
algorithms. The approach for allocating the power within a group is 
to drop the constraint, estimate the quantities of interest and then im- 
pose the constraint via a subsequent normalization. The group-based 
power allocation is computed by 



as,k[i] = Ps,k[i - l]ds. 



k\t 



where 



ks,k 



ds.kli] = ads,k[i] + Zk[i]vk[i\, 
_ a~'^Ps,k[i - l]2;fc[»]ffc[i] 

' l + a--Lv^[{\PsAi~'^]i'k\i]\zk[i\ 



(16) 

(17) 
(18) 



'5.fcW 



~'Ps.k[i - 1] - a~'zl[i\ks,k\i]vi![i\Ps,k[i - !]■ 

(19) 
The normalization as,k[i\ -^ Pa o,s.k[i]/\\&s,k[i\\\ is then made 
to ensure the power constraint. The receive filter is computed by 



Wk[i\ = Pk[i\idk[i\ - p^[{\%'^[{\ip"\i]Pk[i\dk[i] - v) 
where 7"^^ = pf [i]Pfe[i]pj^[i] and 

dk[i] = adk[i] + Zk[i]r[i], 
a-'-Pk[i-l]zk[i]r[i] 



k = 



1 



[i]Pk[i-l\r[i]\zk[i] 



(20) 

(21) 
(22) 



Pfc[i] = a'^Pk[i - 1] - a-\l[i\k[i\r"[i\Pk[i - 1]. (23) 

The proposed scheme employs the algorithm in ( 114) to allocate the 
users in the group and the channel estimation approach of (115b. The 
alternating optimization strategy uses the recursions (I16b -(I23) with 
1 or 2 iterations per time instant i. 



5. SIMULATIONS 

The bit error ratio (BER) performance of the proposed blind joint 
power allocation and interference suppression (BJPAIS) scheme and 
algorithms with group-based constraints (GBC) is assessed. The BJ- 
PAIS scheme and algorithms are compared with blind schemes with- 
out cooperation (BNCIS) 1 16] and with cooperation (CIS) using an 
equal power allocation across the relays (the power allocation in the 
BJPAIS scheme is disabled). A DS-CDMA network with randomly 
generated spreading codes and a processing gain A'^ = 16 is con- 
sidered. The fading channels are generated considering a random 
power delay profile with gains taken from a complex Gaussian vari- 
able with unit variance and mean zero, L = 5 paths spaced by one 
chip, and are normalized for unit power The power constraint pa- 
rameter PA,k is set for each user so that one can control the SNR 
(SNR = PA,k/(T'^) and Ft = Pa + {K - G)PA,k, whereas it fol- 
lows a log-normal distribution for the users with associated standard 
deviation equal to 3 dB. The DF cooperative protocol is adopted and 
all the relays and the destination terminal use linear CCM receivers. 
The first experiment depicted in Fig. [T] shows the BER perfor- 
mance of the proposed BJPAIS scheme and algorithms against the 
BNCIS and BCIS schemes with n^ = 2 relays. The BJPAIS scheme 
is considered with the group-based power constraints (BJPAIS- 
GBC). All techniques employ RLS-type algorithms for estimation 
of the channels, the receive filters and the power allocation for each 
user. The results show that as the group size G is increased the pro- 
posed BJPAIS scheme and algorithms converge to approximately the 
same level of the cooperative training-based JPAIS-MMSE scheme 
reported in [11], which employs G — K for power allocation, and 
has full knowledge of the channel and the noise variance. 



N=16, SNR = 15 dB,K=6 users, P=1500 



^ BNCIS(n =0) 

B BCIS(n^=2) 

■3(6 BJPAIS-GBC(n=2,G=1) 

^ BJPAIS-GBC(n=2,G=3) 

^ BJPAIS-GBC(n=2,G=K) 

JPAIS-GPC[11] 




500 1000 1500 

Number of received symbols 
BER pertormance versus number of symbols. Parameters: p 



Fig. 1 

Ay = Afe = 0.025 (for MMSE and CCM schemes), a 
0.017 and Pfe[j] = 0.01/. 



: 0.998, P. 



s.kW 



The proposed BJPAIS-GBC scheme is then compared with a 
non-cooperative approach (BNCIS) and a cooperative scheme with 
equal power allocation (BCIS) across the relays for n^ = 1, 2 relays. 
The results shown in Fig. |2]illustrate the performance improvement 
achieved by the BJPAIS scheme and algorithms, which significantly 
outperform the BCIS and the BNCIS techniques. As the number 
of relays is increased so is the performance, reflecting the exploita- 
tion of the spatial diversity. In the scenario studied, the proposed 
BJPAIS-GBC with G = 3 can accommodate up to 3 more users as 
compared to the BCIS scheme and double the capacity as compared 
with the BNCIS for the same BER performance, without the need 
for training data. The curves indicate that the GBC for power alloca- 
tion with only a few users is able to attain a performance close to the 
BJPAIS-GBC with G = K users, while requiring a lower complex- 
ity and extra network signalling. A detailed study of the signalling 
requirements will be considered in a future work. 




v:„ 1 ncD .SNR(dB) Number of users ^, . , 

r ig. Z. BbR pertormance versus SNR and number or users tor the optimal 

linear MMSE detectors. Parameters: p = 1, a = 0.998, Ps,fc[i] = O.Oli" 

andPfe[i] =0.01/. 



6. CONCLUSIONS 

This work has proposed the BJPAIS scheme for cooperative DS- 
CDMA networks with multiple relays and the DF protocol. A CCM 
design for the receive filters and the power allocation with group 



constraints has been devised along with a blind channel estimator 
and RLS-type algorithms. The results have shown that the BJPAIS 
scheme achieves significant gains in performance and capacity over 
existing schemes, without requiring training data. Future work will 
consider distributed space-time coding and synchronization. 
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